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ABSTRACT 

Relaxed dark-matter haloes are found to exhibit the same universal density profiles 
regardless of whether they form in hierarchical cosmologies or via spherical collapse. 
Likewise, the shape parameters of haloes formed hierarchically do not seem to depend 
on the epoch in which the last major merger took place. Both findings suggest that 
the density profile of haloes does not depend on their aggregation history. Yet, this 
possibility is apparently at odds with some correlations involving the scale radius Tg 
found in numerical simulations. Here we prove that the scale radius of relaxed, non- 
rotating, spherically symmetric haloes endowed with the universal density profile is 
determined exclusively by the current values of four independent, though correlated, 
quantities: mass, energy and their respective instantaneous accretion rates. Under this 
premise and taking into account the inside-out growth of haloes during the accretion 
phase between major mergers, we build a simple physical model for the evolution of 
Ts along the main branch of halo merger trees that reproduces all the empirical trends 
shown by this parameter in TV-body simulations. This confirms the conclusion that 
the empirical correlations involving rg do not actually imply the dependence of this 
parameter on the halo aggregation history. The present results give strong support to 
the explanation put forward in a recent paper by Manrique et al. (2003) for the origin 
of the halo universal density profile. 
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1 INTRODUCTION 

High-resolution A''-body simulations of hierarchical cosmolo- 
gies with standard cold-dark-matter show that the spheri- 
cally averaged density profile of relaxed dark haloes in the 
present epoch is universal. It is always well fitted by the 
analytical function 



p(r) 



Per a 



rl^{rs -|-r)3-/5 



(1) 



where r is the radius, /3 the central asymptotic loga- 
rithmic slope, with fi xed value equal to 1 accordi n g to 
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2003) or closer to 


1.5 according to other authors (e.g.. 


Moore et al.l Il998l 


iGhigna et alj2000HFukushiee & Makinok200Li:lKlvDin et alJ 


120011). and pc and the scaling density and radius. Thus, 



halos with a given mass have density profiles characterised 
by one single parameter. In fact, the scale radius or the 
concentration c, defined as R/vs in terms of the virial ra- 
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dins R, are found to tightly correlate with halo mass in all 
the cosmologies studied by NFW, a relationship that is be- 
lieved to reflect the fact that less massive haloes typically 
form earlier, when the mean cosmic density is higher (NFW; 
Salvad or-Sole, Solanes & Manrique 1998). iBuUock et alJ 
J200lfl . Wechsler et al. (2002, hereafter WBPKD), and Zhao 
et al. (2003a, hereafter ZMJB) have both extended these xe- 
sults to other redshifts in the flat A CDM cosmology (see 
also lEke. Navarro fc Steinmetd 1200 ll : Zhao et al. 2003b; 
iTasitsiomi et alJl2004l) ' and provided toy-models for the evo- 
lution of rs. 

Within the hierarchical paradigm of structure forma- 
tion, mergers appear to contribute decisively to the evolution 
of haloes. For this reason, it seems natural to try to explain 
the origin of the universal halo density profile f ocusing on the 



role o f mergers (e.g. . Salvador-Sole et al. 1998:ISver fc Whitt 



1998; 



1998; 



(Ra ig. Gonzalez-Cas ado fc Salvador-Sole 

I Subramanian. Cen fc Ostrikerl 
^ ' 



gooo 

a number of m- 



vestigations lAvila- Reese, Firmani & I 


ernanded 


1998 


Nusser fc Shethril998l; iDel PodoIo et al.l 


2OOOI: iKuUl 


200c 


Williams. Babul fc DalcantonI |2004) have demonstrated 



that NFW-like density profiles can be obtained by simple 
spherical collapse, it is not clear why this functionality 



2 Salvador-Sole, Manrique & Solanes 



should be preserved after major mergers since these events 
produce an important loss of memory through the rear- 
rangement and violent relaxation of the system. Yet, as 
shown by numerical experiments, haloes end up with very 
similar density profiles regardless of whether or not they 
have undergone major mergers (|Huss. Jain fc Steinmetzl 
[WQ^) and, provided they have, regardless of the epoch 
when the last of these events took place (iMoore et alJ 
[l999; WBPKD). The situation is even more confusing 
when one takes into account that the above mentioned 
mass-concentration relation and other correlations involving 
rs (or c) recently reported by WBPKD and ZMJB (see 
§ 2J suggest, on the contrary, that the structure of haloes 
depends on their individual aggregation^ history. 

In a recent paper, Manrique et al. (2003, hereafter 
MRSSS) have shown that the sole hypothesis that haloes 
grow inside-out with the typical accretion rate of the cos- 
mology under consideration automatically leads to a den- 
sity profile of the NFW-like form as well as to the mass- 
concentration relation found at 2 = by NFW in nu- 
merical simulations of that cosmology. Moreover, the mass- 
concentration relation predicted for the ACDM universe at 
different redshifts is also in good agreement with the empir- 
ical correlatio ns obtained by Bullock et al. (2001) (MRSSS; 
iHiotelidliooi) . MRSSS argued that the only reasonable ex- 
planation for the good predictions of this pure inside-out 
accretion scheme is that the typical density profile of re- 
laxed haloes is essentially determined by their current mass 
and the boundary conditions set by the matter that is cur- 
rently being accreted. In this manner, relaxed haloes emerg- 
ing from major mergers would have radial profiles indistin- 
guishable from those of haloes with identical global properties 
and boundary conditions but having endured a more gen- 
tle growth. T h is wo uld also explain the fact pointed out 
by iGao et aD i2004) that some aspects of massive ellipti- 
cal galaxies imply that haloes evolve by preserving their 
inner structure as in t he pure inside-out accretion scheme 
jLoeb fc Peebles! I2OO3II while they actually experience im- 
portant rearrangements in major mergers. However, this 
possibility is apparently in contradiction with the correla- 
tions involving r^ (or c) found in numerical simulations. 

In the present paper we examine this issue in detail. 
In §|51 we show that, under some common approximations, 
Ts depends exclusively on the current values of the halo to- 
tal mass, energy and the respective instantaneous accretion 
rates, so there is no room for a dependence on the halo 
aggregation history. From this result, we develop, in § 13 
a simple physical model for the evolution of the scale ra- 
dius along the main branch of halo merger trees that al- 
lows us to show, in § 21 that there is actually no contradic- 
tion between the independence of rs from the halo aggrega- 
tion history and the correlations found in numerical simu- 
lations suggesting the opposite conclusion. Our results are 
summarised and discussed in § |^ Unless otherwise stated, 
the cosmology used throughout this paper is the same as in 
WBPKD and ZMJB studies, that is, a ACDM model with 
{n^,9.K,h,as) = (0.3,0.7,0.7,1). 



^ Throughout the present paper, we use the word "aggregation" 
to refer to halo growth through mergers of any strength and re- 
serve the word "accretion" for growth due to minor mergers only. 



2 GENERAL DEPENDENCE OF THE SCALE 
RADIUS 

We will assume hereafter that relaxed haloes are non- 
rotating and spherically symmetric (see § |S| for a discus- 
sion on the implications of these first-order approximations) . 
Taking the potential origin such that ^(-R) = —GM/R, the 
virial relation adopts the simple scalar form 



dM(r; 



M(r) 



-2E + A-K R'^P 
G 



(2) 



where G is the gravitational constant, E the total energy of 
the halo, P the pressure at the virial radius R and M{r) the 
mass within r, 



M{r) = 



dr 4nf'^ p{r) . 



(3) 



Since the virial radius of a halo at a given cosmic time 
t is a function of its total mass through the relation 



R 



3M 



47rAv(t)p(t) 



1/3 



(4) 



resulting from the identification of haloes as connected re- 
gions with overdensity Av(t) relative to the mean cosmic 
density p{t), it is clear from equation ||5J that p(r) depends 
on the three independent quantities M, E and P. 

Moreover, as we will see in § relaxed haloes grow 
inside-out between major mergers, which allows one to ex- 
press P in terms of the instantaneous mass and energy accre- 
tion rates. Indeed, under these circumstances their structure 
within the instantaneous virial radius at any time t during 
smooth accretion remains frozen. Then the definition of the 
total mass (eq. |3| with r equal to R) and the virial relation 
(eq. HI) lead to 



M{t) - Mo 



fl(t) 



dr 47rr p(r) 



(5) 



2E{t) - 2£o = - 



R(t) 



, 2 . ^ GM{r) 
4Tvr p(r) 



dr 



Ro 



+ A7V [R^{t)P{t)~RlPo] 



(6) 

where R{t), M{t), E{t) and Pit) are the virial radius, mass, 
energy and confining pressure of the evolving halo at t and 
-Ro, Mq, Eq and Po their respective values at an arbitrary 
initial epoch to after the last major merger. By differentiat- 
ing these two equations and taking into account the Jeans 
equation dP/dr — —GM{r)p{r) /r"^ , we obtain the relations 

M = 4-KR'{t)p[R{t)]R 



E = M 



GM(t) 
'~R(tr 



I 3P(t) 1 



(7) 
(8) 



from which one can infer, aside from th e inside-out evolving 
density profile (see MRSSS and lHiotehaii2003l . the pressure 
at the instantaneous virial radius 



Pit) 



M 



6TiK^{t)R 



E GM{t) 



(9) 



where _R is a known function of M, M and t following from 
the differentiation of equation (0J with R and M some func- 
tions of t. 
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Thus, the density profile p{r) of relaxed haloes at t de- 
pends on the four independent quantities: M, E, M and 
E, fixing the right hand-side member of equation Q. Of 
course, these are not the only variables on which p{r) may 
depend since the uniqueness of the solution of equation ^ , 
regarded as an integral equation for p(r), is not warranted 
in general (not even in the case that p(r) can be restricted, 
once M is fixed, to a one-parameter function). Therefore, 
p(r) may in principle depend on other quantities, perhaps 
related to some aspects of the aggregation history of the halo 
as the correlations mentioned in §0 suggest. In Appendix 
1X1 we show however that for density profiles of the NFW- 
like universal form the solution of equation ||5J is unique. 
Consequently, the only parameter required to fix the den- 
sity profile in this case, r^, must depend exclusively on M, 
E, M and E. In other words, all relaxed (spherical, non- 
rotating) haloes with identical current values for these four 
independent quantities necessarily have the same value of r^, 
regardless of their individual aggregation history. 

A first implication of this result is that it allows one to 
understand the origin of the mass-concentration relationship 
shown by haloes at any given epoch without the need to 
presume any dependence of r^ on the time of their formation. 
As we have just shown, rs (or c) depends, at a given t, on M, 
E, M and E, while the latter three quantities are expected to 
correlate, in any given cosmology, with M. To see this latter 
point we must take into account that the kinetic energy 
of protohaloes at some early enough, otherwise arbitrary, 
time ti arises essentially from the Hubble flow at that epoch, 
whereas their potential energy should not be far from that 
of a peak of density contrast 



on the scale 



3M 



1 1/3 



An p{ti 



(10) 



(11) 



with the most probable density profile (|Bardeen et al.ll98(fl . 
In equation (IIOII . 5c is the linearly extrapolated density con- 
trast for spherical collapse at the present time (equal to 
1.69 in an Einstein-de Sitter universe) and D{t) the, also 
cosmology-dependent, linear growth factor. Thus, once the 
arbitrary time ti is fixed, equations (llOt and IIH can be 
used to estimate the (conserved) total energy E of haloes as 
a function of M and t. On the other hand, the mass accre- 
tion rate M of haloes at t is well approximated by (Raig, 
Gonzalez-Casado, & Salvador-Sole 2001) 



ra(M,t) 



AM rps (M, AM, t) d AM , 



(12) 



where rps(M, AM, t) is the usual Extended Press-Schechter 
(EPS) merger rate |Lgx;ov & Cole 1993.) and Am is the frac- 
tional mass increase above which a merger is considered ma- 
jor and, therefore, does not contribute to accretion anymore 
(Salvador-Sole et al. 1998). Finally, by introducing the for- 
mer two functions, E{M,t) and M[M,i), into the funda- 
mental relationship 



one can also estimate the energy accretion rate E in terms 
of M and t. 

We want to stress that this interpretation leads to quan- 
titative predictions of the mass-concentration relation from 
the EPS theory that are, except for rare very massive haloes, 
in good agreement with the results of Ai'-body si mulations 
for different cosmologies and redshifts (MRSSS; iHiotehj 
l2003ft. 



3 A PHYSICAL MODEL FOR THE 
EVOLUTION OF 

We have just seen that the mass-concentration relation can 
be explained without presuming any dependence of the den- 
sity profile of haloes on their mass aggregation history. In § 
|1] it will be shown that the same is true with regards the cor- 
relations reported by WBPKD and ZMJB. But to do this 
we need first to determine how Ts evolves along the main 
branch of halo merger trees. 

In §|Hit was assumed that, during the gentle accretion 
phase between major mergers, relaxed haloes grow inside- 
out, i.e., their inner structure remains frozen. This implies 
that haloes evolving along typical accretion tracks, repre- 
sented by the solution of the differential equation (Raig et 
al. 2001) 



— =r,[M{t),t] 



(14) 



(13) 



dM 
At 

where rs,{M,t) is the mass accretion rate given by equation 
1121 1. maintain an unaltered value of rs or, equivalently, that 
such typical mass accretion tracks coincide with typical rg- 
isopleths. Moreover, according to the result of §|21 these ra- 
isopleths do not depend on the past evolution of haloes but 
just on their current mass. Thus, haloes undergoing a major 
merger will jump from the rs-isopleth they were initially 
following to a new r-s-isopleth according only to their change 
in mass, with the constant value of rs associated with each 
accretion track given by the mass-concentration relation at 
any arbitrary redshift. Therefore, the typical value of rs for 
a halo evolving along the main branch of a merger tree is 
simply that corresponding to the typical rg-isopleth (or mass 
accretion track) that is intersected, at each moment, by the 
Mass Aggregation Track (MAT) followed by the halo. 

At this stage, it is important to remark that the key as- 
sumption in this scheme, namely that haloes grow inside-out 
between major mergers, is fully supported by the results of 
numerical simulations. When a halo evolves through accre- 
tion and major mergers, M necessarily increases but rs may 
not do so; in periods of inside-out growth, rs will be kept 
unaltered. The opposite is not true, of course: the invariance 
of rs does not guarantee that the inner structure of the halo 
remains frozen as the characteristic density pc may vary. 
What should certainly be a secure signature for the inside- 
out growth of haloes along the main branch of merger trees 
is the simultaneous constancy of both the scale radius rs and 
the mass Ms interior to it, as these two quantities fix the two 
degrees of freedom of the NFW-like density profile. In this 
case, the inner structure remains necessarily frozen and the 
mass M grows just because R does. Therefore, the inside- 
out growth of haloes between major mergers must lead to 
values of rs and AIs that are kept simultaneously constant 
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along the main branch of halo merger trees between any two 
sudden increases of M marking major mergers. As Figures 2 
and 3 of ZMJB illustrate, this schematic behaviour is indeed 
found in numerical simulations (see also our Figs. |^ and |S| 
below where we depict the results of Monte Carlo simula- 
tions built under the assumption of halo inside-out growth 
during the accretion phase s). It is no t inco nsistent either 
with the results obtained bv lGao et al.l i2004) from numeri- 
cal simulations of cluster-sized haloes showing that most of 
the matter that ends up, at z = 0, near the centre is not at 
the centre of their main progenitor at z = 4. The reason is 
that the inside-out growth only holds between major merg- 
ers, while the last major merger where current cluster-sized 
haloes were rearranged took place at a redshift much smaller 
than 4 (typically at z ~ 1). 

From Figures 2 and 3 of ZMJB (and Figs. El and |S1 be- 
low) it is also apparent that the phases of inside-out growth 
are only found at relatively low redshifts. There is noth- 
ing strange in this behaviour. The inside-out growth of a 
halo during spherical infall is strictly warranted provided 
only that the orbital period of particles in any given shell 
is much smaller than the characteristic time of mass infall 
(see MRSSS and references therein) or, equivalently, that 
accretion is slow enough. In periods of fast accretion, shell- 
crossing is expected rather to induce the violent relaxation 
of the halo just as major mergers do. As one goes back in 
time, the orbital period of particles in the outermost shell 
becomes of the order of, and even larger than the typical 
infall time of accreted matter in the cosmology considered, 
so the conditions for inside-out growth no longer hold. The 
redshift Za at which the transition between fast and slow 
accretion regimes typically occurs can thus be estimated by 
comparing the orbital period of particles at the virial radius 
of the relaxed halo, 



-7rGAv(t)p(t) 



-1/2 



(15) 



with the characteristic time of mass infall during accre- 
tion, equal to M/ra(M, t). Thus, given the equality p{t) = 
3-ff^(f)r2m(t)/(87rG), where H{t) and 0.ui{t) are respectively 
the Hubble and density parameters at t, and that, in the 
relevant range of cosmic times, f2m(t) and Av(t) can be re- 
sp ectively approximated b y 1 and 200 (for Av(t) a ccording 
to 'Brva n fc NormanllToQi. for which ~ 0.21: iHioteiy 
[2003), the condition for slow accretion adopts the form 



ra[M(f),t] < lQH{t)M{t). 



(16) 



In Figure we show the redshifts Za at which the typical 
mass accretion rates encountered along MATs corresponding 
to different halo masses satisfy the equality rJ[M(t),t\ = 
1.2H{t) M(t), where the factor 1.2 has been chosen to give 
the best agreement with the results of ZMJB (see S 14.21 . 

According to the previous discussion, the typical r^- 
isopleths strictly correspond to the typical evolutionary 
tracks followed by haloes during slow accretion. That is, to 
infer them we must solve equation 1141 taking into account 
only the slow contribution to the accretion rate. Therefore, 
for z ^ Za where accretion is typically slow, the typical r^- 
isopleths will coincide with the solutions of equation 11411 
based on equation II12II . while, for z > Za where accretion is 
typically fast and the contribution of slow accretion to the 
right hand-side member of equation 1141 is null, they will 



14 



13 - 



P 12 h 



11 



10 




Figure 1. Typical mass accretion rates (solid lines) compared 
with the quantity 1.2 H(t) M(t) (dashed lines) found along the 
typical main branch of halo merger trees, here referred to as 
MATs, leading at z = to halo masses equal to 10^^, 10^^, 10^"', 
and 10^^ Mq (from bottom to top) for the same cosmology used 
by ZMJB and WBPKD. The typical MAT corresponding to each 
halo mass is drawn from equation (18) for the appropriate value 
of a according to WBPKD. The epoch during which each pair of 
curves intersects (big dots) gives the characteristic redshift sep- 
arating the fast and slow accretion regimes in the corresponding 
MATs. 



be kept frozen to the values reached at Za (see Fig. I^J. In 
this manner, our model for the typical evolution of rs along 
MATs becomes valid at all redshifts. Note that, compared to 
other prescriptions for the evolution of Vs available in the lit- 
erature, our model stands out for being physically motivated 
and for eliminating the requirement of knowing (and, hence, 
of storing) the individual aggregation histories of haloes in 
order to infer the typical value of their rs at any moment 
since this is estimated from their current mass only'^. 

Given that the inner structure of haloes, determined 
by the values of rs and Ms, is kept unaltered along slow 
mass accretion tracks, the typical rs-isopleths coincide, of 
course, with the typical Ms-isopleths. The constant value of 
Ms associated to an rs-isopleth can be obtained from the 
constant value of this latter parameter and the mass of the 
halo at any arbitrary time along that track from the relation 



Ms = 



(In 2- 1/2) Af 



ln[l + R{M)/n] - [i?(M)/rs]/[l + R{M)/r,] ' ^^^^ 
arising from the NFW universal density profile, with R{M) 



A FORTRAN programme that computes the typi- 
cal rs value for relaxed haloes in CDM cosmologies ac- 
cording to our prescription can be downloaded from 



[http : / /www . am ■ u b . es/cosmo/gravitation.htmj 
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Figure 2. Typical rs (and Ms) isopleths (solid lines) for haloes 
with the same final masses (from bottom to top) as in Figure 
1. The Ta values associated with the different curves have been 
calculated from t he ma ss-concentration relation at z = given 
by^uUock ct a^ 120011) . (The corresponding values are given 
by eq. 1171 .) The curved part of each isopleth corresponds to the 
slow accretion regime, while the horizontal part corresponds to 
fast accretion. Dashed curves show the typical MATs of haloes 
for the same final masses. 



given by equation Q. (Hereafter we assume eq. [1] with 
(3 = 1; were /3 different from unity, eq. |17| should be re- 
placed by the corresponding expression.) This coincidence 
between the rg- and Mg-isopleths implied by our model leads 
to an interesting prediction that can be easily checked. The 
relation between rg and Ma along the different isopleths de- 
fines an Va — Ms relationship that should be typically sat- 
isfied by haloes regardless of their mass and redshift. Such 
a relationship can be obtained from equation 11711 and the 
mass-concentration relation at any arbitrary redshift. There- 
fore, we are led to the conclusion that the r^ — Ma relations 
inferred from the mass-concentration relations at different 
redshifts must overlap. (Note that this conclusion is inde- 
pendent of the form of the ra- or Ms-isopleths and the ex- 
act value of /3.) In Figure |3| we plot the ra — Ma relations 
that result from the empirical mass-concentration relations 
obtained by Zhao et al. (2003b) at different redshifts, as 
well as those calculated, for the same redshif ts and ACDM 
cosmo logy, from the toy- models provided bv lBullock et all 
i200lfl and Eke et al. (2001). In the latter two cases, the 
ra — Ms relations at different redshifts are all very similar, 
indeed, confirming our expectations. There is only a slight 
shift with z which can be attributed to the inherent inaccu- 
racy of the toy-models. However, in the case of Zhao et al., 
the rs — Ms relations show a substantial ^-dependence. This 
discrepancy between the behaviours of the rs — Ms relations 
derived from the mass-concentration relations provided by 
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Figure 3. Typical rs — Ma relations obtained from the mass- 
concentration relations at different redshifts drawn from Zhao et 
al. (2003b) data (top) and from the toy models provided by Eke 
et al. (2001) (center) and Bullo ck et alJ EoOl) (bottom) for the 
same ACDM cosmology as in Zhao et al. (2003b), identical to 
that used throughout the rest of the paper (and in WBKPD and 
ZMJB studies as well) but for the normalisation constant erg, here 
equal to 0.9. 



different authors is not unexpected given the notable de- 
viation of the mass-concentration relations of Zhao et al. 
from those of the other two groups at high z (see Fig. 2 of 
Zhao et al. 2003b). The origin of such a deviation is unclear 
and deserves further investigation. We note however that 
the Zhao et al. mass-concentration relations traced by the 
points drawn from their 25 Mpc size box simulation, 
with softening length equal to 2.5 Kpc, closely follow 
those of the other two groups inferred from simulations with 
similar softening lengths (equal to 1.8 and 0.4 Kpc in 
Bullock et al. and Eke et al., respectively), while the discrep- 
ancy begins to be apparent and becomes definitely marked 
when the Zhao et al. mass-concentration relations are traced 
by the points drawn from their 100 and 300 Mpc size 
box simulations, with softening lengths respectively equal to 
10 and 30 /i"^ Kpc. 
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Figure 4. The c — a(2f) correlation obtained from 2000 Monte Carlo MATs derived in the framework of the EPS formalism using our 
physical model for the evolution of rs for haloes with masses (at 2 = 0) larger than 1.43 X lO^'^ M©. The quantity af plotted in the 
right panel is the normalized cosmic expansion factor, a(zf)/a(0), corresponding to the redshift of formation, Z{, of each halo obtained in 
the way prescribed by WBKPD (see text) from its individual MAT, whereas in the left panel it is the typical value (mean identified by 
circles and median by crosses) for haloes in final mass bins of 200 objects. The straight lines show the best linear fit to the correlation 
obtained from TV-body simulations by WBPKD (solid line) as well as its shift (dashed line) so to fit the correlation obtained by the same 
authors from Monte Carlo simulations using their toy-model for r's that assumes some dependence of this parameter on the individual 
halo MAT. 



4 REINTERPRETING SOME EMPIRICAL 
CORRELATIONS 

4.1 The c vs. a correlation 

WBPKD showed that the MAT^ of individual haloes of mass 
Mo at a redshift zo (as well as the average MAT of haloes 
with Mo at zq) is typically well fitted by the function 



M{z) = Mo exp <^ -a 



a{zo) 
a{z) 



(18) 



with a{z) the cosmic scale factor at the redshift z. Thus, by 
defining the redshift z{ of formation of haloes as the epoch 
in which the logarithmic slope of their respective MAT takes 
some fiducial value S' (~ 2), the previous result implies the 
following relation between Zi and a 



a = —r^S. 

a[zo) 



(19) 



In addition, these authors found that the concentration c of 
haloes at a given redshift correlates with their redshift of 
formation or, equivalently, with the parameter a associated 
to their individual MAT. From this result, WBKPD con- 
cluded that the density profile of haloes is determined by 
their aggregation history. 

There is a different possibility, however. As shown in §|21 
mass and energy conservation, together with the fact that, in 
any given cosmology, the typical accretion rate of haloes is a 
function of Af and t, allows one to relate the concentration c 



^ WBPKD did not distinguish between minor and major merg- 
ers and used the word accretion for what is here referred to as 
aggregation. 



of haloes at some given redshift to their mass. This implies, 
in turn, that c also correlates with any arbitrary function of 
M, in particular with a{M) giving the typical (mean or me- 
dian) value of a for haloes of mass M. Since the individual 
values of a obviously scatter around their typical one, the 
c — a(M) correlation arising from the mass-concentration re- 
lation will induce a secondary correlation between c and the 
a value corresponding to each MAT or, equivalently, with 
its associated a{z{) value. In other words, the c — a{zf) corre- 
lation reported by WBPKD might simply be a consequence 
of the well-known mass-concentration relation and, just like 
it, might not imply any dependence of rs on the aggregation 
history of haloes. 

To check the validity of this alternative explanation we 
ran 2000 Monte Carlo MATs for presen t -day haloes with 
masses according to the lSheth fc TormenI (1 19991) mass func- 
tion and satisfying the same restriction, M > 1.43 x 10^^ 
Mq, as in WBPKD and followed, in each, the evolution of 
Ts predicted by the model developed in § |3 that does not 
presume any dependence of this parameter on the halo ag- 
gregation history. These Monte Carlo MATs were built from 
the conditional probability that haloes of a given mass at 
a given epoch arise from s maller objects at ea rlier epochs 
given by the EPS theory jLacev fc Colelll993l) . The time 
steps adopted were small enough to ensure the inclusion 
of only one merger with AM above the ma s s reso lution of 
the simulations according to iLacev fc Col j il993l) . In this 
way, to be sure that the inferred MATs really correspond 
to the main branch of merger trees, we only had to impose 
that the relative mass increase in any time step was smaller 
than unity. The values of Ta corresponding to the different 
isopleths in our physical model for the evolution of this pa- 
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rameter were drawn from the mass-concentration relation at 
z = Q of iBullock erai] (1200 ll) . 

In Figure]^ we show the c — a{zf) correlation at z = 
obtained from these Monte Carlo simulations. To compen- 
sate for the fact that our model for Vs does not include ran- 
dom deviations of E, M and E from their typical values 
that depe nd ex c lusive ly on Af and t, we have added, ac- 
cording to Ijind ||2pOCf), Gaussian deviates in log(c) with a 
dispersion equal to 0.09. The main trend of this correla- 
tion reproduces satisfactorily that found in WBPKD's N- 
body simulation (see their Fig. 6). The only small differ- 
ences concern the slightly larger scatter and the small shift 
(a little more marked in the medians than in the means) 
in the Monte Carlo correlation relative to the A'^-body one. 
This behaviour is identical to that shown by the correlation 
obtained by WBKPD in the Monte Carlo simulation they 
performed, using a slightly different implementation of the 
EPS theory, to check the validity of their toy-model for 
that includes an explicit dependence of this parameter on 
the individual halo MAT. These authors found, indeed, a 
scatter systematically larger (by 25%) in the Monte Carlo 
correlation than in the A''-body counterpart and the corre- 
lation itself shifted exactly by the same amount as found in 
the present work. As mentioned by these authors, the small 
deviations of the Monte Carlo predictions with respect to 
the Af-body results are likely caused by the well-known in- 
accuracies of the E PS formalism (and its pra ctical imple- 
mentation; see e.g. ISomerviUe fc Kolatli [19991 and Raig et 
al. 2001). In any case, what is important here is that our 
physical model for Vs is capable of recovering the correlation 
found by WBPKD at the same level as their toy-model de- 
spite the fact that, unlike their model, ours does not presume 
any dependence of Vs on the halo aggregation history. 



4.2 The Ms vs. correlation 

ZMJB recently found that haloes show a tight correlation 
between the Vs and Ais values scaled to those of the main pro- 
genitor at a redshift ztp where the MAT shows a characteris- 
tic change of slope** , which they interpreted as an indication 
of the passage from fast to slow aggregation regime. More- 
over, they showed that the log- log slope of this scaled Vs — Ms 
correlation changes depending on whether the halo redshift 
is larger or smaller than the respective Ztp value. This find- 
ing was interpreted as evidence that the inner structure of 
haloes depends on their aggregation history. Still, ZMJB did 
not offer any explanation neither for the origin of this cor- 
relation nor for its dependence on the halo redshift relative 

to Ztp. 

In fact, the frequency of major mergers cannot be held 
responsible for the change of slope in the MAT of a halo 
at Ztp because, in the ACDM cosmology, major mergers 
are, for the present-day halo masses considered by ZMJB 
(~ 10" Mq), even more frequent in the slow than in the 

* The discontinuity in the MAT slope coincides with a maximum 
or "turning-point" in the function log(?;c) -I- 0.25 Iog[-ff (z)], where 
Vc is the circular velocity of the halo along its MAT, which pro- 
vides a well-defined practical procedure to determine ztp. Note 
also that, like WBPKD, ZMJB used the word accretion for what 
is here referred to as aggregation. 



fast aggregation regimes defined according to the heuristic 
redshift ztp. What could explain, instead, the characteristic 
change of slope in MATs observed at that redshift is the 
transition from fast to slow accretion at Za mentioned in 
§ Besides, the correlation between rg and Ms predicted 
by the model developed in that section is expected to split 
into two tight correlations similar to those found by ZMJB 
when the values of the two quantities are scaled to those of 
the main progenitor at z^. Note that the — Ms relation 
shown in Figure|3|is very nearly linear in log-log. Therefore, 
the fact that major mergers and fast accretion tend to in- 
crease with decreasing redshift the values of rs and Ms as 
haloes evolve along MATs plus the slight z-dependence of 
the slope of the log-log rs — Ms correlation derived from the 
Zhao et al. (2003b) mass-concentration relations (see Fig. 
1^ will translate into a small bending (on top of the slight 
intrinsic one; see Appendix B) of the log-log rs — Ms correla- 
tion shown by haloes along MATs. Consequently, by scaling 
the values of rs and Ms to those of the main progenitor at 
Ztp (or 2a) along each individual MAT, such a bent corre- 
lation will split into two essentially linear relations on both 
sides of that redshift, with that corresponding to the fast 
regime being slightly steeper as found by ZMJB. Note how- 
ever that, according to our model, since the values of rs and 
Ms do not depend on the particular MAT the haloes have 
followed, the scaled rs — Ms correlations reported by ZMJB 
should not rely on their scaling by the specific values of rs 
and Ms corresponding to the main progenitor found at ztp 
(or 2a) along each particular MAT. Using as reference the 
object found at Ztp (or 2a) along the typical MAT for haloes 
with the mass of that being studied should yield essentially 
the same results. 

The validity of our interpretation of the ZMJB findings 
can be checked, like in the case of the WBPKD correlation, 
by means of Monte Carlo simulations. Unfortunately, this is 
not so immediate in the present case. As explained in §|21 our 
model for the evolution of rs (and Ms) is not fully consistent 
with the mass-concentration relations found by Zhao et al. 
at high redshifts. However, we can still check it on the scaled 
rs — Ms correlation arising fr om the mass-concentration rela- 
tions of lBullock et al.1 J200lll or Eke et al. (2001) that do not 
suffer from this problem. As, in the latter case, the general 
rs — Ms correlation does not depend on z, we do not expect 
its 2-induced bending along MATs. But, as shown in the 
Appendix^ this general correlation still has some intrinsic, 
less marked, convexity, which will also make it split, when 
properly scaled, into two linear correlations similar to those 
found by ZMJB, with the reference halo in the individual or 
typical MATs playing exactly the same role. 

With this aim, we performed Monte Carlo MATs as ex- 
plained in § 14.11 but now following also the evolution of AIs . 
In Figures 13 and |S| we show a few of these realisations for fi- 
nal halo masses similar to those appearing in Figures 2 and 
3 of ZMJB. To better mimic the ZMJB results, the mass 
increases produced along the MATs were accumulated be- 
tween each couple of consecutive output times chosen iden- 
tical to those in the ZMJB A'^-body simulations. The ztp 
values, marked on these tracks by arrows, were obtained 
following the prescription given by ZMJB, whereas the 2a 
values, marked by big dots, were calculated using the same 
criterion as in Figure Despite the fact that our rs and 
AIs tracks do not include the effects of the random devia- 
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Figure 5. Monte Carlo realisations (upper curves) of MATs and the corresponding Ms tracks (lower curves) derived from the values 
of Tb plotted in Figure |B] for the same cosmology and similar final halo masses as in ZMJB. The location of the redshift Za (big dots) 
marking the typical frontier between fast and slow accretion regimes is compared to the redshift Ztp (indicated by arrows) marking the 
characteristic knee in the mass aggregation tracks indicated by ZMJB. To compare with Figure 2 of ZMJB. 



tions of E, M and E from their respective typical values, 
they show a remarkable similarity to those depicted in Fig- 
ures 2 and 3 of ZMJB. There is good agreement not only 
in the general shape of the MATs and the associated and 
Ms tracks (see, in particular, the effects of inside-out growth 
during slow accretion), but also in the typical location of the 
Zfp values. From these figures we can see that every zt-p is 
close to the corresponding Za marking the change of accre- 
tion regime along each MAT, which confirms our guess on 
the origin of the change in slope of MATs. It is interesting to 
note that the procedure followed to determine Ztp tends to 
locate this characteristic redshift just after a major merger 
or an interval of fast accretion. This introduces a systematic 
discontinuity in the slope of the MATs at that redshift that 
explains the two-branch analytical fitting formula proposed 
by Zhao et al. (2003b) for the average MAT of haloes. How- 
ever, when Ztp is not privileged as in WBPKD, the typical 
MAT of haloes of mass M shows no such discontinuity. 



In Figure |7| we plot the scaled — Ms correlations 
corresponding to the fast and slow accretion regimes, as in 
Figure 7 of ZMJB, drawn from objects along the Monte 
Carlo MATs of Figures|^and|S|at the same redshifts as in the 
ZMJB figure. To mimic the effects of the random deviations 
of E, M and E in A''-body haloes, we have added once again 
Gaussian deviates in log(c) with a dispersion equal to 0.08 
and 0.12 in the slow and fast accretion regimes, respectively. 
(Unlike WBPKD, who impose no restriction on their halo 
sample, ZMJB chose relatively isolated final haloes. Hence, 
the dispersion in c of these haloes is expected to be slightly 
smaller at z < Za than in the case of WBPKD. On the 
contrary, the fact that haloes should be typically less relaxed 
at z > Za should cause a higher dispersion in c there.) In the 
top panels, we show the results obtained when the scaling 
reference is the halo located at Za along the typical MAT 
for the mass of each object, while in the bottom panels we 
show those obtained by using, instead, the halo at ztp along 
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Figure 6. Evolving virial radii R (upper curves) and scale radii Tb (lower curves) predicted by our physical model for the evolution of 
this parameter corresponding to the Monte Carlo MATs plotted in Figure|F] To compare with Figure 3 of ZMJB. 



each individual MAT just like in ZMJB. As expected, both 
variants show identical trends: tight linear correlations in 
log-log with similar slopes and scatters, strongly resembling 
those in the original ZMJB scaled correlations. Specifically, 
the slopes in the fast and slow accretion regimes are 0.76 ± 
0.03 and 0.65±0.03, for the typical-MAT variant, and 0.79± 
0.02 and 0.64 ± 0.03, for the particular-MAT one. Since in 
the former case, the correlations obtained do not involve, 
by construction, the MAT of individual haloes, we conclude 
that the change in slope of these scaled correlations at ztp 
(or 2a) does not imply the dependence of Vs (and Ms) on 
the individual halo aggregation history. The same conclusion 
should therefore hold for the correlations obtained by ZMJB 
(with slopes in the fast and slow regimes respectively equal 
to 0.64 and 0.48). 



5 SUMMARY AND DISCUSSION 

In the present paper, we have shown that the scale radius 
Ts of relaxed, non-rotating, spherically symmetric haloes en- 



dowed with the universal density profile a la NFW depends 
exclusively on the current values of their mass and energy 
and the instantaneous accretion rates of these two quan- 
tities. There is therefore no room for any additional de- 
pendence on the particulars of their aggregation histories. 
Using a simple, physically motivated model for the evolu- 
tion of ra along the main branch of halo merger trees, we 
have demonstrated that our finding that rg does not de- 
pend on the halo past evolution is fully compatible with 
the well-known mass-concentration relation and the correla- 
tions recently reported by WBPKD and ZMJB suggesting 
the opposite conclusion. We note however that the correla- 
tion found by ZMJB has an origin slightly different from that 
predicted within our scheme since it appears to be related to 
the particular ^-dependence of the mass-concentration rela- 
tion found by these authors compared to that obtained by 
other groups. The present results support the explanation 
recently provided by MRSSS for the origin of the universal 
halo density profile. 

Apart from the correlations here investigated, N- 
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Figure 7. Scaled rg — Ms correlations drawn from objects along the Monte Carlo MATs shown in Figures |F] and [5] at the same output 
times of the ZMJB simulations for the slow (left) and fast (right) accretion regimes. The values r^^o and M^^q used to scale the and Ms 
variables and to separate both aggregation regimes correspond to the halo found at along the typical MAT of each object (top) and 
to the halo found at ztp along the real individual MAT of each object (bottom). Different symbols are used for objects along different 
MAT realisations. The straight lines show the slopes of the similar, though with a different origin (see text), empirical correlations found 
by Zhao et al. (2003b), shifted by the same amount as in their Figure 7. 



body simulations show that the density of the local 
environment also correlates with the frequency of ma- 
jor mergers occurring along the halo aggregation histo- 
ries (the higher the density, the h igher this frequency; 
iGottlober. Klvpin fc Kravtsovl 1200 j) and with the value 
of Ts (the higher the density, the smaller r„: Ijind l2000l : 
iBuUock et aP 2001). Hence, some correlation may exist be- 
tween Ts and the frequency of major mergers along the main 
branch of individual halo merger trees. As is well-known, in 
its original form the EPS formalism does not contain in- 
formation on the spatial distribution of haloes, so it is not 
possible to infer fr om our model (unle ss modified along the 
lines proposed bv lMo fc Whitelll996l or ISheth etaPbOOlD 
any correlation between Va and local density from our model. 
However, given that the environment of a halo should af- 
fect its instantaneous accretion rate (the higher the density. 



the larger the typical accretion rate) and that, as we have 
demonstrated, the value of Vs depends on this latter quan- 
tity, a correlation between Vs and local density is also ex- 
pected to hold within our scheme. For the same reason, the 
fact that haloes located in high (low) density environments 
may have had a higher (lower) chance of merging with other 
haloes along their past history is not inconsistent with the 
present results either. What these results do not support 
is the existence of a direct causal relation between the fre- 
quency of major mergers endured by a halo during its growth 
and the final value of its scale radius. 

The conclusions of the present study rely on some hy- 
potheses that deserve some discussion. The least controver- 
sial of these is that of the inside-out growth of haloes be- 
tween major mergers; as shown in §13 it is fully supported 
by the results of numerical simulations. In contrast, the as- 
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sumption that haloes are spherically symmetric and non- 
rotating is actually not very realistic. Had we not adopted 
these simplifying approximations, the angular momentum, 
the asphericity of the system (via the inertia tensor) and 
the anisotropy of the local velocity tensor would appear in 
the full vectorial form of the virial relation, invalidating the 
proof in Appendix^ We thus cannot discard the possibility 
that, under more realistic assumptions, the inner structure 
of haloes can preserve some memory of their aggregation 
history. We are not concerned here, of course, about the 
memory embodied in the final value of any extensive prop- 
erty, such as the total mass, energy or angular momentum, 
that are integrated, and hence fully convolved, over the en- 
tire evolutionary history of the system. The kind of memory 
that is relevant for the present discussion is rather that as- 
sociated with the possible imprint on the inner structure 
of haloes of specific events, such as the formation of the 
halo (whatever its definition), the transition from fast to 
slow aggregation regimes or the frequency of major merg- 
ers, having taken place at some point in the past. Certainly, 
the elongation and the anisotropy of the velocity tensor of 
haloes indicate that the episodes of violent relaxation do 
not proceed until completion. But this does not necessar- 
ily mean, of course, that the loss of memory produced in 
those relaxation processes is so small that precise informa- 
tion on those specific events can be recovered from the final 
properties of the spherically averaged density profile. The 
results obtained in the present paper rather suggest, along 
the lines of thos e drawn from A'^-body expe r iments men- 
tione d in § Q ( Huss. Jain fc Steinmet j Il999l : iMoore et alJ 
I1999I : WBPKD), that the density profile of haloes is largely 
insensitive to the details of their aggregation history. In any 
case, they show that, from the correlations found in numer- 
ical simulations, it cannot be concluded that rs depends on 
the halo aggregation history since, at least under the as- 
sumption of spherical symmetry and null-rotation, they are 
well explained and quantitatively recovered without presum- 
ing such a dependence. 
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APPENDIX A: UNIQUENESS OF THE NEW 
SOLUTION OF THE VIRIAL RELATION 

Consider two steady spherically symmetric, non-rotating, 
self-gravitating systems with identical radius R and total 
mass M. This obviously warrants that their respective den- 
sity profiles p(r) and p'{r) coincide at least at one radius 
in the open interval {Q,R). We next prove that if the two 
systems also have identical total energy E and are subject 
to the same external pressure P then the radius at which 
both profiles coincide cannot be unique. 

Suppose that the density profiles p(r) and p'(r) (as- 
sumed with realistic central logarithmic slopes greater than 
—2 to ensure finite central potentials) do intersect each other 
at a single radius within (0,i?). The scalar virial relation 
allows one to write 



AM{r) 



(Ai) 



where M(r) and M' (r) are the mass profiles associated with 
p(r) and p'{r), respectively. By integrating by parts both 
sides of this equality, we obtain 



dr 



dr 



M'^{r) 



(A2) 



Rewriting equation \K2l in the form 



dr 



K^ir) M'^{r) 



dr [M(r) - M'(r)] 



Af(r) + M'(r) 







(A3) 



and integrating by parts the second member leads to 

rR 

dr [p(r) - p'(r)] J(r) = , (A4) 



where we have introduced the definition 

M(c) + M'(e) 



J(r)^ 



(A5) 



Equation IIA4II can be recast in the form 



drr^ [p('') ^ p'('")] J{''') 



drr^ [/9'(r) — /9(r)l J(r) , 



(A6) 



with the integrands on both sides of equation (IA6II being 
two strictly positive (or negative) functions in the respective 
open intervals (0,rc) and {re,R). Since J(r) is a strictly 
outwards increasing positive function, J(re) is, at the same 
time, an upper bound for r € [0, re] and a lower bound for 
r £ [re,-R]. Thus, we have the following inequalities 



Jo 

f 

Jr. 



drr^ [p(r) — p'(r)] J(r) < 
drr^ [p(r) — p'(r)] J(re) 



drr^ [p (j) — p(r)] J(r) > 



(A7) 



drr^ \p(r) - p(r)] J(rc) , 
which, given equation 1A6II . lead to 

2 

J(re) / drr [p(r) - p (r)] > 

J(re) / drr^ [p'(r) - p(r)] 
or, equivalently, to 

drr^p{r)> I Arr^p'(r). 



(A8) 



(A9) 



(AlO) 



As inequality (lAf 01 contradicts the fact that the two den- 
sity profiles must yield the same total mass M, the initial 
postulate that they coincide at one single radius r^ can be 
rejected. 

Two profiles of the NFW-like form, equation (f), yield- 
ing the same mass M, are either identical or intersect at 
one single intermediate radius r^ in the interval (0,7?). In- 
deed, if they are not identical, they have different values of 
rs and of the characteristic density pc, implying that they 
have different slopes at all radii. This ensures that the func- 
tion p(r) — p'(r) has no extreme in the interval (0, i?) and, 
consequently, that it only vanishes once. Thus, according to 
the previous proof, we conclude that, if at a given epoch 
two (spherical, non-rotating) relaxed haloes have identical 
values of R, M, E and P, their density profiles of the NFW 
universal form must necessarily coincide, too. 



APPENDIX B: THE Ms VS. SCALED 
CORRELATION 

Let us neglect, for simplicity, the deviates of E, M, and E 
of individual haloes from their typical values dependent on 
M and t. As explained in § |3 the rs — Ma relation can be 
inferred from the mass-concentration relation at any given 
redshift and equation 11711 arising from the universal density 
profile of hal oes. The mass-co ncentration relation sd, z — 
given by uBuUock et al.l i200ll) for the same cosmology as 
in WBPKD and ZMJB leads to the nearly hnear r^ - Ms 
relation in log-log units shown in Figure IBll Should this 
relation be strictly linear, the rs and Ms values of haloes 
would satisfy the relation 



log Ms = Qo + 3ai log rs 



(Bl) 



Since the values rs,o and Ms,o of any arbitrary halo even- 
tually used to scale the rs and Ms values of the different 
objects would also satisfy the relation IIBIL we are led to 
the fact that haloes also show a perfect linear relation (with 
no scatter) of the form 

log(Ms/Ms,o) = 3ai log(rs/rs,o) . (B2) 

However, the real rs — Ms relation shown in Figure lHTI shows 
a small bending so that haloes rather satisfy the relation 

\2 



log Pis = Qo + 3ai log rs + 3a2 (log rs) 



(B3) 



being ai > 0, a2 < and \a2 log(rs,max/rs,min)| < ai, 
where rs,max and rs,min are the upper and lower bounds of 
rs in the sample. Given that the values rs,o and Ms,o of 
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any reference halo used to scale and Ms also satisfy the 
relation IIB3L we are led to the scaled relation 

log(Ms/Ms,o) = 3ai log(rs/rs,o) 

+ 3a2 [(log rsf - (log rs,of] , (B4) 

which can be rewritten in the form 

log(Ms/Ms,o) = 3ai (^1 + 2 ^ logr^.o) log(rs/rs,o) 

+ 3a2[log(rs/r,,o)]'. (B5) 

Since ]q2 log(rs/rs,o)| ^ tJi, the relation <B5ll is well ap- 
proximated by the linear relation 



log{Ms/Ms,o) = 3ai log(rs/rs,o) 



1 + ^ 

Ql 



21ogrs,o + log 



(B6) 



where rs^m is some intermediate value between r^fi and 
fs.max or ra,min depending on whether is larger or smaller 
than rs,o, respectively. Thus, the slope depends, in this case, 
on the specific value of rg,o used to scale r^. Consequently, 
by taking rs,o equal to the scale radius of some halo along 
the typical or individual MAT of each object, the slope will 
slightly change, yielding some scatter in the scaled rela- 
tion. Note however that this scatter should be very small 
as \a2 log(rs,m/T's,o)/Qi I is always much smaller than unity. 
On the other hand, since rs,ni/rs,o is then respectively either 
smaller or larger than unity, given the signs of coefficients 
ai and 02, the slope of the scaled correlation will be sys- 
tematically larger, though by a small amount, for Vs ^ rs,o 
than in the case rs > Vs.o- 

Finally, we do not expect any noticeable difference be- 
tween the Za, and Ztp variants of the scaled correlations. As 
mentioned in f0| the values of along individual MATs are 
of the same order as Ztp, so their respective values of rs,o 
should be similar. Certainly, the redshifts correspond to 
objects along typical MATs, but individual MATs are not 
very different from the latter. Thus, both variants of the 
scaled correlations should show essentially the same slope 
and scatter. 



QO 
O 




-7 -6 
3 log(ryMpc) 

Figure Bl. Typical Tb — Ms relation obtained from the same 
haloes as in Figures|5]and|n](same symbols as in Figure|7} without 
adding any scatter in c and prior to the scaling of and Alg 
and the subsequent splitting in two relations on both sides of 
the reference halo. To better appreciate the small bending of this 
relation we show its log-log linear fit (full line) shifted vertically 
by 0.07. 



